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NOMENCLATURE 


c’^(S) n times continuously differentiable on S 
p pressure 

U component velocity in x-direction 

V component velocity in y-direction 

A[ ) + aM ]/3^[ ]/3y^ 

AA[ ] 3H ]/ax^ + 2 + g**[ j/gy** 

ax"' ay'^ 

as boundary of domain S 
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SUMMARY 


For more than the last 20 years there has been a concerted effort to 
solve the stationary Navier-Stokes equations; however, this has only been suc- 
cessful for a few special cases of primarily academic interest. An alternative 
approach has been to solve the equations numerically, and then compare the 
results with experiment. On occasion, such comparisons are in good agreement. 
However, such results are of dubious value since one has no a-priori way of 
knowing the relevance of such results until they are explicitly compared 
against experiment. Therefore, it would seem reasonable to conclude that the 
present approaches to solving the Navier-Stokes equations are of limited 
value. Accordingly, it is the purpose of this paper to show that there does, 
indeed, exist an equivalent representation of the problem that has significant 
potential in solving such problems. This is due to the fact that this equiva- 
lent representation of the problem consists of a sequence of Fredholm Integral 
Equations of the second kind, and the solving of this type of problem is very 
well developed. In addition, for the problem in this form, there is an excel- 
lent chance to also determine explicit error estimates, since one would now be 
dealing with bounded linear operators, rather than unbounded. 


INTRODUCTION 


Ideally, one would ]ike to obtain a numerical solution to the Stationary 
Navier-Stokes problem that is within some prescribed degree of accuracy of the 
true solution. Unfortunately, it is not possible to accomplish this objective 
by just applying some existing numerical technique directly to the Stationary 
Navier-Stokes equations, since error analysis for nonlinear equations is, for 
all practical purposes, nonexistent. However, since it is possible to replace 
the Stationary Navier-Stokes problem by an equivalent sequence of linear 
partial differential equations (ref. 1), providing some rather general condi- 
tions can be met, it is at least theoretically possible to obtain error esti- 
mates due to the vast wealth of knowledge known about linear equations. How- 
ever, even though attaining such error estimates is within the realm of 
possibility, one should not be misled into thinking that such a task is a 
small undertaking, for this is most certainly not the case. However, the 


chances of someone necompUshing this task nvo significantly greater If the 
basic seqvjentlnl problem Is transformed into an equivalent form that is more 
suitable to error analysis. Accordingly, in this paper it will be proven that 
the. said sequential problem can be transformed into a form that has signifi- 
cantly greater changes of yielding eKplleit error estimates. 


GENERAL DEVELOPMENT 


The Navior-Stokes equations are given by 

uu^ + vUy *f p^ - V An + fj(x,y) * 0 

V 


uvj^ + vvy + p^ - V Av + fgCx.y) » 0 


in S 


0) 


(2) 


with boundary conditions 

u(OS) « -b^OS) , v(3S) - bj , 

whore S is a two dimensional Green's domain with surface 9S, 

fj(x,y) f C^(S) , and f,(x,y) p C^S) , 

or equivalently, 

vAAt!' + Ai; - li* A?' + f- -f, «0 in S 

>A JV) iy 4 .^ 

V„(9S) ■ b, , ij* (DS) « b, 

•* y 

Let 

vAAv'' + >1 AV„ “ Ai|' + fj^ - fg 
y X X y ly ‘-J, 

be denoted by PC','). Hence (2) can bo expressed by 

PC';) - 0 . I'^Os) - bj , VyCas) « b. 

where P can be interpreted as a mapping from c‘'(S) into C^fS). As demon- 
strated in reference 1, solving (2) is equivalent to solving the sequence of 
equations 


P(»|,j ) + P' (i|i ) (v'^ 1 , I/' ) « 0 


3x 




9S 


1 » 


!!d 

9y 


(3) 


as 


for r. « 0,1 ,2 , . , . , 


2 



providing certain conditions are satisfied (ref, 1), (i.e., primarily that 
tfti the initial guess is reasonably good)» for further insight into the sig- 
nificant latitude on i|!q> the render should consult reference 2. 


3 * vAAt ] + A[ + Ayj, 




Equivalent to equation (3) Is the problem 


3$ 

P(«n> + P' - 0 . -X 


n 


0 

8s ■ 


0 , n « 0| 1 ,2, 


» W 


by letting 




n+i 


Now consider (A) under a slight change in boundary conditions, in particular: 

3?' 


P(i|< ) + P'(i|i ):jj - 0 , ij» 
'^m'^ '^m ra 


3S 


0 , 


3n 


3S 


0 , m « 0,1,2, . . 


(5) 


If is a solution to (5), it directly follows that is also a solution 


of (A)» and vice versa (i.e., equation (A) li 


m 


3S 


constant, without loss of 


generality assume constant « 0). Therefore, it is sufficient to focus our 
attention on equation (5), which in detail is given by 


With 


iilj„(3S) » 0 , (3S) “ 0 , m “ 0,1,2, . . . 


Hence, 


where G is the Green's function of the biharraonic Equation [3] 
AAi}) “0 in S, with M3S) ■= 0 , i{)^(3S) = 0 


Lemma: 


(6) - (5) 


Proof. D 


3 


( 6 ) 


for con'^cnionee and elnrityi oquaLion (6) will be oxpresBed in the form; 
where 

> ■ V Jg M 1„ + A«n^l ly - '"n.j,! 1» - S 

and 

b 

Lemma; If the hypothesis of theorem 1 of reference 1 is satisfied, then there 
exists a solution of equation (6) for all m ond 

n 

y,* - am E “ 'I’o 

iyho i«o 


where V'*'' is the unique of equation (2). 

Theorem; Under the condition of the above lemma, is a solution to the 
Fredholm Integral Equation of the second kind 


♦*(x' 


I 


y') + J 'l•J(x,y)KJ,,(x' ,y' ,x,y)dx dy 


f('i„,)(x',y') 


(7) 


where 


K„,(x’ ,y' ,x,y) » -AG^ (x' ,y* ,x,y)iim^^(x,y) + + 2 G^.^ 


+ 'ImijjyC^xx “ Sy^ 


"xy 

Proof; Since ^ solution of (6) under the hypothesis, it is suffi- 

cient to show that a solution to (6) is alos a solution of (7). for clarity, 
denote (6) in the fonn 


" f(vo) 


or equivalently, 


- X 

V 


f g(, 0 ,* + 


AV’o Jl'jj - 'i’o 


^ Atj)ds « f(Vo) 




A 


By Greon’s Lemma 

"J* G'I'o dx dy « J A fo p * dx dy 

s y s ' y 

v)x'^ * " ['^(^'* “y)]' 

" I 4\)f 

“/ G M ' 0 V'J d& « +J (g AV- ) ij,* ds 
*^s X y •'s'' x'v 


‘3S 


I 


+ I G ds 


■/ 

^ r« 


+J Gt|u Aig ds 
s K y 


-J (g a*; 0 ) y 

s'^ y'x 

■/ i^\)j 


Thorefoi*e (6) reduces to 


<!-=*+-/ 1 ’*(a(gv -'5 ) + (g ) - (g ) Ids - f(te) 

'S ' y'x ' x'y ' x'y ' y'X’' 

^\\ - 4 \\ - «\)J 

« d(Gjji*-'o + Gi^ o ) + G A^o + G Ai|»o - lIg A ^ + Gs^ „ ) - G M‘o ~ G A<; ^ 

\ *' y xy/ ' X xy N ' M xy/ y x 

** ^(px'^o ” ^y'*'o (^y "" ^x 

Kvo ) “ (M'o ) + (g,^^o J 

' y' ' y'xx ' y'yy 

- °xxx»»^ + “x't'O + tixI'D^ + C^yy»0 + «vv) 
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'{VoJ • ■ iv\\y 


(■ 




) <■ 

kX/v ' 


^yy’^0^ " 


xy'y 


”®xxv'*'o ” ^v'^0 “ ^^xv'*’u ** ^y'^0 " ^vvy** 0 '"^®yy**o 

XKy Ojj y Oj,yy xy y y;.y 


^(“x*o ■ V't^ ■ ^°xy('' 0y ■ ♦•xx)'*' ^’''xy(°’‘’‘ " *'’'i 


+ ♦. 40x - ♦, 4Cy 

y * 


.’. i^(pjjV o - ^Gy A’t'o “ Gx ^*' 0 ) 

y / 

“ ^[^xy(^0 “ ’^Oxx) ■*■ '•'«xy(°'''' " ’*’ ^'®y ^ 


Therefore, (6) reduces to 

'j* V / '^*['^0 ^*^x ~ ’^'0 ^^y ^(^Ky(^o “ ’^%x)’*^'^®xy(^^^” 

s y ^ ^ ' 

Hence, (6) *> (7) . 


Gjjx” Gyyjl ds » f (i, o) 


CONCLUSION 


Therefore, under the conditions as cited, the sequence of solutions 
generated by (7) converge to the solution of (1) in the following sense: 

t it ■ („ - '•■, 

i»c 


E 


i|,« - ^ 


Therefore, 


iji* » ijjo + ]C 


0 ■ ^ ^i 

0 
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Henec the original problem, equation (I), which falls within the frame- 
work of nonlinear operator theory, an area Chat little is known, has been 
replaced by equation (5), which falls within the framework of bounded linear 
operator theory, an area for which there exists a vast wealth of information. 

Of course, it could be argued that it is theoretically possible that the 
sequential representation of the solution converges so r.lowly that the results 
are of questionable value. However, from all Indications this will not be a 
problem, for as pointed out in reference 1, there are parameters at one's 
disposal in the method that can bo adjusted to speed up the rate of conver- 
gence; in fact it was demonstrated in reference 2 that by judicious selection 
of the variable parameters available they were able to get numerically adequate 
convergence with just a few iterations. 
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